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Abstract. In this note, we determine the ground state energy of the trans- 
lation invariant PauH-Fierz model to subleading order O{o^) with respect to 
powers of the finestructure constant a, and prove rigorous error bounds of 
order O(a^). A main objective of our argument is its brevity. 



1. Introduction 

We study the translation invariant Pauli-Fierz model describing a spinless elec- 
tron interacting with the quantized electromagnetic radiation field. We present a 
very short and simple method to determine the subleading terms of the ground state 
energy up to order 0{a^), where < a <C 1 denotes the finestructure constant, 
and to rigorously bound the error by 0{a'^). 

A well-known difficulty connected to this problem arises from the fact that the 
ground state energy is not an isolated eigenvalue of the Hamiltonian, and that the 
form factor in the interaction term of the Hamiltonian contains a critical frequency 
space singularity (the infrared problem of Quantum Electrodynamics (QED)). One 
of the most striking consequences is that the ground state energy does not exist as 
a convergent power series in a. It is inaccessible to ordinary perturbation theory, 
and can only be written as a convergent expansion in powers of a of the form 
^nbn{a) a"^ , with a-dependent coefficients bn{a) that diverge in the limit a ^ 
at a sub-power rate, i.e., |6„(a)| = o{a^^) for any (5 > 0, [T1[31[S]. 

It is possible to determine the ground state energy and the renormalized electron 
mass to any arbitrary precision in powers of a, with rigorous error bounds, by use 
of sophisticated rigorous renormahzation group methods, [ll[3l[5]. However, these 
algorithms are highly complex, and explicitly computing the ground state energy 
to any subleading order in powers of a is a voluminous task. 

Therefore, it remains a desirable goal to devise alternative methods that produce 
such results more directly for non-trivial, but reasonably low orders in powers of a 
with rigorous error bounds, and with much less effort. 

The present paper is the first in a number of works investigating an approach to 
such problems based on perturbations around the true ground state of the interact- 
ing system. A key ingredient in our argument is a new estimate on the expected 
photon number in the ground state derived from [H| whose proof involves a bound 
on the renormalized electron mass uniform in the infrared cutoff, [2 [S] . 

Estimates on the ground state energy play an important role, for instance, in 
binding problems, e.g., the determination of the hydrogen binding energy. The 
leading term of the ground state energy was demonstrated to be of order O(a^) (up 
to normal ordering) in [8] , and explicitly determined, with a rigorous error bound of 
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order 0{a^). A similar result was subsequently obtained for the spin ^ case in In 
[7j, the ground state energy for the Bogoliubov-transformed, translation- invariant 
Nelson model is explicitly determined up to O(a^), with an error bound of order 
0{a^) (uniformly in the ultraviolet cutoff). While for a fixed ultraviolet cutoff, 
this result is fully analogous to ours, we present here a new and particularly short 
method, and our sharp error bound of order O(a^) is necessary for a companion 
paper, where we address the hydrogen ground state energy. 

2. Definition of the model 

We study a non-relativistic electron interacting with the quantized electromag- 
netic field in Coulomb gauge. The Hilbert space accounting for the pure states 
of the electron is given by L^(R'^), where we neglect its spin. The Fock space of 
transverse photon states in the Coulomb gauge is given by 

where the n-photon space = (S)" <Xi C^) is the symmetric tensor prod- 

uct of n copies of one-photon Hilbert spaces L'^{M.^^) (g) C^. The factor accounts 
for the two independent transversal polarizations of the photon. On ^, we intro- 
duce creation- and annihilation operators a^(fc), ax{k) satisfying the distributional 
commutation relations 

[ax{k),al,ik')] = 5x,x'5{k~k') , [a{ik) , a{,{k')] = 0, 

where a\ denotes either a\ or a\ . There exists a unique unit ray flf E "S, the Fock 
vacuum, which satisfies a\{k) fi/ = for all A; S and A e { + , — }. 

The Hilbert space of states of the system consisting of both the electron and the 
radiation field is given by 

n := l2(R3) (g, ^. 

We use units such that h ^ c ^ 1, and where the mass of the electron equals ni — ^. 
The electron charge is then given by e = ^/a, where the finestructure constant a 
will here be considered as a small parameter. 
The Hamiltonian of the system is given by 

H lei Hf+ : {iVx If - VaA{x))^ : (1) 

where ; ( • • • ) : denotes normal ordering (corresponding to the subtraction of a 
normal ordering constant proportional to a), 

Hf = I dk\k\al{k)ax{k) 

is the Hamiltonian for the free photon field, and the operator 

^«^E_I^«<i^i.M.(.>-^«..,. 

couples the electron to the quantized vector potential. The polarization vectors 
s\{k), A £ {+)—}) are unit vectors orthogonal to one another and to fc g M^, in 
accordance with the Coulomb gauge condition, divA = 0. 
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The function k implements a fixed ultraviolet cutoff on the wavenumbers k. 
Via scaling, we may assume k to be compactly supported in < 2}, monotone, 
and to satisfy k — I for |A:| < f . For convenience, we shall write 

A{x) = A-{x) + A+{x), 

where 

is the part of A{x) containing the annihilation operators, and = {A~{x))* . 

The system is translation invariant, and H commutes with the operator of total 
momentum 

Ptot = iVx (E) If + lei ® Pf, 

where iV^ and 

Pf = J2 [ dkkal{k)ax{k) (2) 

denote the electron and the photon momentum operators, respectively. It there- 
fore suffices to consider the restriction of H to the fiber Hilbert space Tip = ^ 
corresponding to the value P e R'^ of the conserved total momentum, given by 

H{P) = Hf+ : {P- Pf -V^A{0)f : . (3) 

Henceforth, we will write A± = A±(0). From [2 [5], it is known that inf (j{H) = 
inf a{H(0)), and we will in the sequel only consider P = 0. 

3. Statement of the Main Results 
On 5^, we define the positive bilinear form 

{v,w). := {v,{Hf + Pf)w), (4) 

and its associated semi-norm — {v,v)l/^. 



Theorem 3.1 (Ground state energy of H{0)). We have 

11*2 11^ 

^a^2\\A-<^2r - - 4||$i||?) + 0{a^) (5) 



inf spec (i/(0)) = -a^ II $2 II* 



wh 



ere 



$2 
$3 

$1 



= -iHf+Pf)-'A+-A+nf, (6) 
- -{Hf + Pf)-'PfA+^2, (7) 
= -{Hf + P])-'PfA-^2. (8) 



Remark 3.1. An easy computation shows that 

C > II $j lU > > (9) 

for a constant C < oo, and j — I, 2, 3. 

Let g 5^ be the minimizer of H{0), normalized by (4" , fi/) = 1. Taking the 
( • , • )*-orthonormal projections of along the vectors j = 1, 2, 3, and denoting 
the component in the ( • , • )*-orthogonal complement of their span by R, we get 

* = il/ + 2 ?7i $1 + 772 a $2 + 2 ?73 $3 + i? (10) 
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where for j — 1,2,3 

<i>,)^ = ||$,||^^., (11) 

= , = = (12) 

The coefficients rjj remain to be determined. 

Theorem 3.2 (Ground state of H{0)). Let * be the ground state (T^ of H{0), 
normalized by (^',51/) = 1, [1, 5J. Then, 

^ ^ + a^2 + 2ai ^1 + 2a^ ^3 + R, (13) 

where 

R := i? + 2 (771 - 1) $1 + (772 - 1) a $2 + 2 (773 - 1) $3 . (14) 

The coefficients rjj satisfy \ 771,3 — 1 1 < ca, and 1 772 — 1 | < ca^ . Moreover, 

\\R\\ , ||i?|| < ca , and \\R\\^ , \\ R\\^ < ca^ . (15) 

A key ingredient in our proof is the following estimate on the expected photon 
number. 

Proposition 3.1. The expected photon number in the ground state is bounded by 
(^', iV/*) < Ca2(^', (16) 
for a positive constant C < 00 independent of a, where 

fdkal{k)ax{k) (17) 



N 



f - 



is the photon number operator. 

Proof. For cr > 0, let Ha{P) denote the fiber Hamiltonian regularized by an infrared 
cutoff implemented by replacing the function k by a function k^- with Kg- — k 
on [cr, 00), Kct(O) = 0, and Kg- monotonically increasing on [0, cr]. Then, Ea-{P) := 
inf spec {H„{P)) is a simple eigenvalue with eigenvector \E'<t(-P) S S^, [HIS]. If P = 0, 
one has VpE„{P = 0) = 0, [1, 5j. In formula (6.11) of [6], it is shown that 

ax(k)^a{Q) = (/) + (//) (18) 

where from (6.12) in [F follows that 

II (7)11 < C{k)\VpE,{0)\ = 0, (19) 

and that 

if the electron spin is zero. Thus, it follows immediately from (6.19) in [B] that 

for spin zero, where mren^a- is the renormalized electron mass for P ~ 0, [HIS], 
defined by 

1 _i ^ (Vp*,(o), (i/,(o) - £;,(o)) Vp*,(o)) 



(vE'.(0),vE'.(0)) 



(21) 
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As proved in [Tl[5], 1 < mren,a < 1 + ca uniformly in a > 0. Correspondingly, we 
obtain 

(^',iV/*) = lim /dfc|laA(fc)*.(0)|l2 < Ca^ll^f 
<^\" J 

as claimed, where = s — limcr\o ^ct(O) (see [1]). □ 

4. Proof of Theorems 13.11 and 13.21 

To prove Theorems 13.11 and 13. 2|, we derive the following upper and lower bounds 
on inf spec {H{Q)). 

4.1. The upper bound. We define the trial function 

^triai ._ _|_ 2a5$i + a$2 + 2a5$3. (22) 
Then, the variational upper bound on the ground state energy 

^ ^trial ^ ^(-Q)^trial ^ _q,2||^^||2 (23) 

+ a3(2||A-$2||' - 4|1<J>3|1J - 4|l$i||2) + Oia"), 
is obtained from a straightforward computation. 

4.2. The lower bound. Substituting the expression pO|) for ^, and exploiting 
{■ , ■ )^ -orthogonality, straightforward calculations explained in Section [5] yield 

(*,i/(0)«') > -a2||$2||' + a'h2p(2M-$2||' - 4||$i||2 - 4||$3||') 
+ a^\m~l\^\\<^2\\l 

+ 4amm~V2\^\\<^i\\l + \m-V2\^\\^3\\l) 

-ca^l + \ru\' + \V2\' + \V3\')- (24) 

We recall from ([9]) that C > || $j ||* > C^^ > for some finite C independent of 
a, and j = 1, 2, 3. We use 

\V2 - Vj\^ > 1\V,^M' - \V2-lf and < 2 - 1 + 2 (25) 

to show that the three last lines in (p4)) are bounded below by 

iTyi-lpa^Bi + |772 -lpa^S2 + |772-lpa^B3 - ca^ (26) 
where the constants Bj > are defined by 

B2 = ||$2||' - 4a(|l$i|l2 + ||$3||2) _ ca^ 
Si = 2||$i||2 - ca 

Bi = 2||$3||'-ca. (27) 
Clearly, (|26p > — ca^, and the minimizing triple (771, ?72j satisfies 

h7i.3-l|<ca , and |?72-l|<ca^ (28) 
Therefore, is bounded from below by 

- «'||'i>2||^ + o?{2\\A-^2f - m^Wl - miWl) + \\\R\\l - ca\ (29) 



Combined with we find || i? ||* < ca\ This proves Theorems O and [321 □ 
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5. Proof of Inequality ([24ll 

Clearly, 

i/(0)vl/) = {yS^,{Hf + Pf)^) (30) 
+ 4:VaRe{-^ , Pf A- (31) 

+ 2ai?e(^', A" A"*) (32) 

+ 2a(^', A+A- . (33) 

We estimate the quadratic form of H{0) on the true ground state ^E". 

5.1. Some auxiliary estimates. Since 

M-^ll < cWHj^PW, (34) 

one has 

M-VII ,\\Pf^\\ < c II VII, (35) 
for all -0 e 3^ in the intersection of the domains of and Pf. Let for brevity 

$ := Qf + 2 r/i $i + ry2 a $2 + 2 r/3 $3 (36) 
so that 5" = 4* + i?. We observe that from ([9]), ([35|) . and the Schwarz inequality, 

\\A-^f, ||P/^ir < ca^\r,2\^ + camml' + . (37) 

In addition, we recall that A^Pj — PfA^ by the Coulomb gauge. 

5.2. The term dSO]). First of all, the term gives 

{^,iHf + Pf)^) = a'\\^2\\l+4am\<f,\\l + II $3 11^) + ||i?||2 (38) 
due to the orthogonality relations pT|) . IT^ . 

5.3. The term ([3T|). We expand the term according to * = * + i?. Clearly, 

+ 8a^Re-n2V3{A+^2, Pf<^>3)- (39) 

Moreover, we have 

4V^Re{R, PfA- R) > -i||i^||^ (40) 

o 

- 

Here, we used || A^ R\\ < c\\ R \ \ combined with ([57]) . and applied the Schwarz 
inequality in the form |^_B| < ^A^ + for any < 5 < 00. 

Furthermore, we find 

8V^Re{^,PfA-R) > -^\\R\\l -ca^irjl + rjl) (41) 
by similar arguments, using 

{<f2,PfA-R) = ($3, i?), = = (flf^PfA-R) (42) 
(see (ini)), combined with ([ST]) . 
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< \\Nf R\\ < 



NJ * II - II A^/ * I 



5.4. The term (|32]). We have 

2aRe{^ , A- A- ^) = 2aRe{A+ ^ , A^ ^) . 

Clearly, 

||A+^f = [A-,A+]||Vf + ||A-^f , (43) 

where < [A^, A+] < c, for all ijj E ^. Using the Schwarz inequality and similar 
arguments as above, one finds 

2ai?e(^', ^- > 2a^ Rer]2{A+ A+flf , ^2) 

-l\\R\\l^ca'\\Rr 

-ca^lmf + + \v^f), (44) 
where we used that if ?/' = ^f, $1, $3, or R, then by ((TTl), (HH), 

{A+A+nf,ij) = ($2, V'), = 0. (45) 
At this point, we apply Proposition 13 . 1 1 in 

' < ca. (46) 

Here, we use the fact that R has a vanishing projection on the Fock vacuum, (H^), 
and II ^ || < ca, which one easily verifies. 

5.5. The term (|33p . This term is estimated by 

2a\\A-^f > 2a3||A-<i>2|p - ca^lml' + Ivsf) - J||i?||^ (47) 

5.6. Collecting all estimates. Combining the 0(a"^)-terms of ([55)) with (|39p . we 
find 

4a^|?7ini$i||J + &a''Re7nm{^i, 

= -Aa^\Tj2\^\\^i\\l + Aa''\vi-V2\^\\'^i\\l, (48) 

since from ([8|), 

($1, P/A-$2) = ~ 11*1 II*. (49) 

Likewise, 

4a3|773|||$3||2 + 8a^Re7j2V3{PfA+^2,^3) 

= -4a3|ry2|'||$3||' + 4a3|,72-%ni$3||L (50) 

since from ([7|), 

(P/A+<I>2, <J>3) - -||<!>3||^ (51) 

Combining the 0(Q!^)-term in ([55]) with the first term on the r.h.s. of 

a'h2ni*2||' + 2a^Rerj2{A+A+nf,<i>2) 

= -«'||<I'2||^ + a2|,y2-l|2||<I,2||2, (52) 

since from ([B]), 

{A+A+nj,<f2) = -||$2||^ (53) 

Moreover, we use half of ||i?||^ in ([55]) to compensate all four terms — 1|| i?||^ ap- 
pearing in the above estimates, and are left with the term i|| i? ||^ in 
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All remaining terms are bounded below by — ca'' ( 1 + + I772P + |?73p )■ 
Collecting all bounds, we arrive at ([24|) . □ 
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